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ABSTRACT 

We identify 13 sets of multiply-lensed galaxies around MACS J0717.5+3745 (z = 0.546), outlining 
a very large tangential critical curve of major axis ~ 2.8' , filling the field of HST/ACS. The equivalent 
circular Einstein radius is e — 55 ± 3" (at an estimated source redshift of z s ~ 2.5), corresponding to 
r e ~ 350 ± 20 kpc at the cluster redshift, nearly three times greater than that of A1689 (r e ~ 140 kpc 
for z s = 2.5). The mass enclosed by this critical curve is very large, 7.4+0.5 x 1O 14 M and only weakly 
model dependent, with a relatively shallow mass profile within r < 250 kpc, reflecting the unrelaxed 
appearance of this cluster. This shallow profile generates a much higher level of magnification than 
the well known relaxed lensing clusters of higher concentration, so that the area of sky exceeding a 
magnification of > 10 x, is ~ 3.5D' for sources with z ~ 8, making MACS J0717. 5+3745 a compelling 
target for accessing faint objects at high redshift. We calculate that only one such cluster, with e > 
55", is predicted within ~ 10 7 Universes with z > 0.55, corresponding to a virial mass > 3 x 10 15 M , 
for the standard AC DM (WMAP5 parameters with 2<j uncertainties). 

Subject headings: gravitational lensing , galaxies: clusters: individual: MACS J0717. 5+3745 , dark 
matter 



1. INTRODUCTION 

Improvements in gravitational-lensing modeling to- 
gether with the fabulous image quality of the Advanced 
Camera have led to the identification of increasing num- 
bers of multiply-lensed images and the discovery of some 
surprisingly large lenses (e.g., Kneib et al., 1996, Broad- 
hurst et al. 2005, Bradac et al., 2008, Halkola et al. 2008, 
Liesenborgs et al., 2008, Limousin et al., 2008, Zitrin et 
al. 2009, Zitrin & Broadhurst 2009). 

The largest bound structures to form hierarchically 
are likely to have collapsed recently and hence should 
be found relatively locally, at low redshifts. The larger 
volume available with increasing distance means that in 
practice we cannot expect to reside next to the most mas- 
sive cluster. The Coma cluster, z = 0.023, is the most 
massive "local" cluster with a reliable weak lensing based 
mass of M vir ~ 0.65 x 1O 15 M (Gavazzi et al., 2009), 
similar to earlier dynamical estimates ~ 0.8 x 1O 15 M 
(The & White, 1986, Geller et al., 1999). Other more 
distant clusters are known from lensing to be more mas- 
sive, such as A2218 at z = 0.16 (Kneib et al., 1996) and 
A1689 at z = 0.18, with 1.6 + 0.2 x 1O 15 M (Broadhurst 
et al., 2005, Umetsu & Broadhurst, 2008), superseded by 
the most distant Abell cluster, A370 at z = 0.37 and the 
most luminous X-ray cluster RXJ1347, z = 0.45, with 
masses reliably determined from weak lensing distortion 
and magnification of ~ 2 x 1O 15 M (Broadhurst et al. 
2008). 

Strong lensing is not seen for low redshift clusters 
(z < 0.15) because the central projected mass densities 
do not exceed the critical density required to generate a 
sizable Einstein radius, which scales inversely with lens 
distance, diverging at low redshift. Similarly, at high 
redshift as the lens distance approaches that of the dis- 
tant sources, thus the critical density is too great for 
strong lensing. Geometrically, lensing is optimized at in- 
termediate redshifts, where for a given mass the critical 



density for lensing is minimal, but this is partly offset by 
the late hierarchical growth of high-mass systems. This 
trade-off results in estimates of the amplitude of strong 
lensing to favor the redshift range z = 0.2 — 0.4, for 
NFW-like mass profiles (Broadhurst & Barkana, 2008, 
Oguri & Blandford, 2009). Seemingly, all rich and X- 
ray luminous clusters at modest redshift are found to 
have many multiply-lensed images when examined with 
sufficient depth and resolution, implying that the mass 
profiles of clusters are in general sufficiently peaked that 
the critical density for lensing is exceeded. 

The Einstein radii of the most massive clusters seem to 
be larger than predicted by the AC DM model (Broad- 
hurst & Barkana, 2008, Puchwein & Hilbert, 2009), based 
on the "Millenium" simulation (Springel et al., 2005). 
This discrepancy is empirically supported by the surpris- 
ingly concentrated mass profiles measured for such clus- 
ters, when combining the inner strong lensing with the 
outer weak lensing signal (Gavazzi et al., 2003, Broad- 
hurst et al., 2005 & 2008, Limousin et al., 2008, Don- 
narumma et al., 2009, Oguri et al., 2009, Umetsu et al., 
2009, Zitrin et al., 2009). 

These largest lensing clusters have proven to be excel- 
lent targets for accessing the faint early Universe due to 
their large magnification consistently providing the high- 
est redshift galaxies (Ebbels et al., 1996, Franx et al., 
1997, Frye & Broadhurst, 1998, Bouwens et al., 2004, 
Kneib et al., 2004, Bradley et al., 2008, Zheng et al., 
2009, see also Broadhurst, Taylor & Peacock 1995). The 
larger sizes of lensing images has provided increased spa- 
tial detail and the large magnifications permits obser- 
vations of increased spectral resolution, leading to the 
discovery that metal enriched material is typically out- 
flowing from galaxies at z > 4 (Franx et al. 1997, Frye 
& Broadhurst, 1998, Frye, Broadhurst & Benitez, 2002). 

With the goal of discovering high redshift galaxies 
and to better define the mass profiles of galaxy clus- 
ters in general, we have combined data for a sample of 
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well studied clusters. In this process we have uncovered 
the extreme lensing properties of MACS J0717. 5+3745 
(z=0.546), a cluster originally identified in the highly 
complete sample of the most X-ray luminous clusters in 
the Universe (Ebeling et al., 2004, 2007). This cluster is 
thought to be amongst the most massive clusters known, 
forming part of an intricate dark matter (DM) filamen- 
tary structure of ~ 4 Mpc (Ebeling et al., 2004, Ma et al., 
2008), whose complex X-ray emission indicates ongoing 
merging (Ma et al., 2009) accompanied by the most pow- 
erful known radio halo (van Weeren et al., 2009, Bonafede 
et al., 2009). 

In §2 we describe the observations; the lensing analysis 
is described in §3; our results are presented in §4 and in §5 
we discuss and conclude them. Throughout the paper we 
adopt the standard cosmology (£l m o = 0.3, Q\o = 0.7, 
h = 0.7). Accordingly, one arcsecond corresponds to 
6.42 kpc/hjQ at the redshift of this cluster. The refer- 
ence center of our analysis is fixed near the center of the 
ACS frame at: RA = 07:17:31.65, Dec = +37:45:03.12 
(J2000.0). 

2. OBSERVATIONS 

The very X-ray luminous cluster MACS J0717.5+3745, 
which is the denser north-western region of the large- 
scale filament found by Ebeling et al. (2004), was imaged 
in April 2004 and in October 2005, 2006, with the Wide 
Field Channel (WFC) of the ACS installed on HST. In- 
tegration times of ~ 4500s were obtained through each 
of the F555W and the F814W filters in the April 2004 
run. We retrieved these images from the Hubble Legacy 
Archive, along with a subsequent exposure in the F606W 
band (January 2005, integration time of 1980s) to form 
a composite high-resolution 3-color image. Several obvi- 
ous close pairs of multiply-lensed galaxies and giant arcs 
are immediately visible throughout the full frame, with 
which we begin our modeling process described below. 

3. LENSING ANALYSIS 

We apply our well tested approach to lens modeling, 
which we have applied successfully to A1689, C10024, 
and MACS Jl 149.5+2223, uncovering large numbers of 
multiply-lensed images in several clusters imaged with 
HST/ACS (Broadhurst et al., 2005, Zitrin et al., 2009, 
Zitrin & Broadhurst 2009). The full details of this ap- 
proach can be found in these papers. Briefly, the basic 
assumption adopted is that mass approximately traces 
light, so that the photometry of the red cluster member 
galaxies is the starting point for our model. 

Cluster member galaxies are identified as lying close 
to the cluster sequence by the photometry provided in 
the Hubble Legacy Archive. We approximate the large 
scale distribution of matter by assigning a power-law 
mass profile to each galaxy, the sum of which is then 
smoothed. The degree of smoothing and the index of the 
power-law are the main free parameters. A worthwhile 
improvement in fitting the location of the lensed images 
is generally found by expanding to first order the gravita- 
tional potential of the smooth component, equivalent to 
a coherent shear describing the overall matter ellipticity, 
where the direction of the shear and its amplitude are 
free, allowing for some flexibility in the relation between 
the distribution of DM and the distribution of galaxies, 
which cannot be expected to trace each other in detail. 



The total deflection field dr(0), consists of the galaxy 
component, d ga i(9), scaled by a factor K ga i, the clus- 
ter DM component ^m^), scaled by (1-K ga i), and the 
external shear component a ex (6): 

a T {0) = K gal a gal (6) + (l-K gal )a DM (0)+a ex (0), (1) 

where the deflection field at position m due to the ex- 
ternal shear, d ex (6 m ) = (a ex , x ,a ex , y ), is given by: 

&ex,x (0 m ) = M cos(20 7 )Aa; m + | 7 | sin(2</> 7 )A?/ m , (2) 

®-ex,y \Vm 

) = |7| sin(2^ 7 )Ax m - |7| cos(2(/> 7 )A?/ m , (3) 

where (Ax m ,A?/ m ) is the displacement vector of the 

position # m with respect to a fiducial reference position, 
which we take as the lower-left pixel position (1,1), and 
<p 1 is the position angle of the spin-2 external gravita- 
tional shear measured ant i- clockwise from the x-axis. 

We lens candidate galaxies back to the source plane 
using the derived deflection field, and then relens this 
source plane to predict the detailed appearance and 
location of additional counter images, which are then 
searched for in the data (see Figures 1 and 2). The fit is 
assessed by the RMS uncertainty in the image plane: 

RMS? magea = ^((x;-x0 2 + (^-^) 2 ) / N image8 , (4) 

i 

where x i and y i are the locations given by the model, 
and Xi and yi are the real images location, and the sum 
is over all N images images. 

Importantly, this image-plane minimization does not 
suffer from the well known bias involved with source 
plane minimization which biasses solutions towards high 
magnification and hence correspondingly shallower pro- 
files. The model is successively refined as additional 
sets of multiple images are identified and incorporated 
to improve the fit (Zitrin et al., 2009). In addi- 
tion, we find that a bright foreground elliptical galaxy 
(RA=07:17:37.16, DEC=+37:44:22.54) is important to 
include in the model as it locally affects lensed images in 
the eastern to central part of the ACS frame, including 
their relative lensing distance. 

In the above process we uncovered and used 34 
multiply-lensed images, corresponding to 13 lensed back- 
ground galaxies, to constrain the mass distribution and 
profile of this cluster. A lensed pair of red drop-out ob- 
jects which are not seen in the F555W band at an es- 
timated redshift of z ~ 4 (system number 5, see Figure 
1), helps us pin down accurately the normalization of the 
deflection field by adopting the lensing distance for this 
redshift which is nearly independent of cosmology (see 
also Broadhurst et al., 2005, Zitrin et al., 2009). The 
image-plane RMS of our final model is very good, 2.2 " 
per image, whereas for comparison, in A1689 an RMS 
of 3.2 " was achieved by Broadhurst et al. (2005) and 
for C10024 an RMS of 2.5 " was achieved by Zitrin et al. 
(2009). 
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Fig. 1. The cluster MACS J0717.5+3745 (z = 0.546) imaged 
with Hubble/ACS F555W, F606W, and F814W bands. The 34 
multiply lensed images identified by our model are numbered here. 
The white curve overlaid shows the tangential critical curve cor- 
responding to the distance of system 1 at an estimated redshift 
of z ~ 2.5, and which passes through several close pairs of lensed 
images in this system. The larger critical curve overlaid in red cor- 
responds to the larger source distance for the red dropout galaxy 
number 5, at the estimated photometric redshift of z ~ 4. This 
large tangential critical curve encloses a very large lensed region 
equivalent to ~ 400 kpc in radius at the redshift of the cluster, 
z = 0.546. 





Fig. 3. — Radial surface mass density profile, «(r), in units of 
the critical surface density, i.e. k(t) = E(r)/E cr ^, derived for the 
range of radius covered by all sets of multiple images shown in 
Figure 2. The profile is shallow within ~ 250 kpc, at a level close 
to the critical density for system 1 (black line). The vertical dashed 
line at ~55" is the Einstein radius at the redshift of system 1. This 
profile was measured circularly while centered on the central pixel 
of the 2D mass distribution shown in Figure 4. 





Fig. 4. — Zoomed-in 2D surface mass distribution («), in units 
of the critical density. Contours are shown in linear units, derived 
from the mass model constrained using 34 multiply-lensed images 
seen in Figure 2. The axes are in ACS pixels. Note that the 
central mass distribution of is rather flat reflecting the unrelaxed 
appearance of this cluster. 



Fig. 2. — A reproduction of system 1 by our model. The observed 
images are shown in the top row are compared with our model 
generated images below. Each model image is generated using 
as input the pixels of image 1.1 (except for the model image of 
1.1 which is generated by relensing the observed image 1.2) and 
delensing these pixels back to the source plane and then relensing 
the source plane to generate the counter images. It is clear that our 
model is successful in demonstrating the multiply-lensed relation 
between the four observed images. 



4. RESULTS 

The derived surface mass distribution is relatively shal- 
low (see Figure 3) in accordance with the unrelaxed ap- 
pearance of this cluster, which is known to be in the 
process of merging, with disturbed X-ray emission, hot 
shocked regions (Ma et al., 2008, 2009), and powerful ra- 
dio halo emission (van Weeren et al., 2009, Bonafede et 
al., 2009). We find that the elongated central distribu- 
tion of galaxies is followed by a very extended tangential 
critical curve, enclosing a large critically lensed region of 
~ 2.63D' with an equivalent Einstein radius of ~ 55±3.0" 
(for an estimated source redshift of z ~ 2.5). This corre- 
sponds to a physical scale of 350 ± 20 kpc/hm, substan- 
tially larger than in any other known cluster. This very 
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Fig. 5. — The number of clusters per universe with masses ex- 
ceeding M, lying at redshift z > 0.55, for which the predicted 
Einstein radius is > 55 ± 3 /; . The shaded area reflects the vari- 
ance generated by incorporating into the the ACDM model the 2a 
uncertainty in the cosmological parameters h, n, and ag, about the 
WMAP5 values, and also the la uncertainty in the measured Be- 
This number peaks near M~4x 1O 15 M , with N c i ~ 2.7 • 10 -7 . 

large radius adds to the already uncomfortable discrep- 
ancy between the large Einstein radii observed for mas- 
sive clusters and the predictions based on the standard 
ACDM cosmology (Broadhurst & Barkana, 2008, Sadeh 
& Rephaeli, 2008, Puchwein & Hilbert, 2009) for which 
such large Einstein radii can only be contemplated with 
mass distributions which are highly prolate and aligned 
along the line of sight (Corless & King, 2007, Oguri & 
Blandford, 2009). Instead, here this cluster is evidently 
elongated across the line of sight, traced by the distri- 
bution of member galaxies, and by the extended X-ray 
emission and the tangential curves (see Figure 1). 

Naturally, such a comparison is biased by the asym- 
metry and shape of the critical curve, which cannot be 
perfectly compared with DM simulations. However, this 
difference is small, as we find that a circularly-symmetric 
lens of0E = 55", centered on the center of Figure 4, con- 
tains 6.8 x 1O 14 M , 92% of the mass contained within 
the tangential critical curve (for z s = 2.5). 

The mass enclosed within the tangential critical curve, 
for z s = 2.5, is very high 7.4 ± 0.5 x 10 14 M o . The 
mass enclosed within the larger critical curve for the 
multiply-lensed dropout galaxy, number 5, at an esti- 
mated source redshift of z ~ 4 is correspondingly larger, 
~ 1.0 x 1O 15 M (marked in red, Figure 1). Clearly, the 
total mass associated with this cluster cannot be smaller 
than this and is likely to be several times larger, and may 
soon be reliably estimated from weak lensing. We stress 
that the scale of the tangential critical curve hardly de- 
pends on the mass profile because it is set by the critical 
density for lensing, which depends only on fundamental 
constants and the distances involved, so that different 
profiles explored in the modeling procedure yield very 
similar critical curves, and hence similarly large Einstein 
radii and an accurate enclosed mass. 

5. DISCUSSION AND CONCLUSIONS 

Theoretically, we can make an approximate estimate 
of the probability of obtaining such a cluster in the con- 



text of the tightly proscribed ACDM model. Using the 
probability distribution function (PDF) of halo forma- 
tion times, calculated within the framework of the ex- 
tended Press & Schechter (EPS) formalism (Lacey & 
Cole 1993), we derived a corresponding PDF of halo 
concentration parameters (Sadeh & Rephaeli, 2008) by 
adopting the formation redshift-concentration scaling, 
deduced by Wechsler et al. 2002, from N-body simu- 
lations. 

Solving the equation governing the relation between 
the concentration parameter and 0e assuming an NFW 
profile (e.g. Broadhurst & Barkana 2008), leads to a 
PDF of Einstein radii. Here we adopt the mean source 
redshift z s = 2.5 as described above. The product of the 
Press & Schechter mass function abundance of clusters 
within the relevant mass and redshift ranges with the cu- 
mulative Einstein radius probabilities multiplied by the 
volume available at z > 0.55 provides the number of clus- 
ters expected above this redshift for standard ACDM 
where we allow for 2a uncertainty in the WMAP5 val- 
ues of the parameters [n, h : as]). Extremely low numbers 
of clusters are predicted peaking at M ~ 4 x 1O 15 M 
with a maximum probability of only 2.7 x 10 _T objects 
per Universe (see Figure 5), or approximately half this 
value when allowing for the area of sky covered by the 
ROSAT all-sky survey. This comparison with ideal DM 
halos is made somewhat uncertain by the observed asym- 
metry of the central mass distribution. However as noted 
in section 4, there is only an 8% difference between the 
mass within the critical area and the mass within the 
circular equivalent critical area derived from our model, 
whereas the discrepancy with theory is orders of magni- 
tude. Clearly it is important to obtain the total mass of 
the cluster via weak lensing which we are now estimating 
but it is reasonable to suppose this lies close to the most 
probable value derived crudely here, given that the mass 
interior to the critical curve is already 25% of this value. 

MACS J0717. 5+3745 is the second strong-lensing ana- 
lyzed cluster (along the recent critically convergent lens 
MACS J1149.5+2223, Zitrin & Broadhurst 2009) out of 
several other very luminous X-ray clusters uncovered by 
Ebeling et al. (2007) in the ROSAT all-sky survey. Both 
these clusters have unconcentrated central galaxy distri- 
butions reflecting their unrelaxed state and for which we 
have shown that the central mass distribution is likewise 
very unconcentrated. A relatively shallow mass profile 
boosts the gravitational lens magnification and we cal- 
culate that the total area of sky exceeding a magnifica- 
tion, [i > 10, is ~ 3.5 square arcminutes, corresponding 
to the current high redshift limit of z ~ 8, which is far 
higher than the equivalent area calculated for other mas- 
sive clusters (Broadhurst et al., 2005, Zitrin et al., 2009). 
Such unrelaxed and massive clusters opens a potentially 
new regime of highly magnifying lenses for accessing the 
faint distant Universe. 
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